r; t=0.01 9:51:17

. version 15.1
r; t=0.00 9:51:17

. clear   all
r; t=0.00 9:51:17

. macro   drop _all
r; t=0.00 9:51:17

. set     linesize 80
r; t=0.00 9:51:17

. set             more off
r; t=0.00 9:51:17

. set             varabbrev off
r; t=0.00 9:51:17

. set     rmsg on
r; t=0.00 9:51:17

. 
. // #0. Document this do-file and set up packages and directories
.         local today : di %tdCY/N/D daily("$S_DATE", "DMY")
r; t=0.00 9:51:17

.         tokenize `today',parse("/")
r; t=0.00 9:51:17

.         local year `1'
r; t=0.00 9:51:17

.         local month `3'
r; t=0.00 9:51:17

.         local day `5'
r; t=0.00 9:51:17

. 
.         local who "Ben Cheung"
r; t=0.00 9:51:17

.         local pgm "20250823-brfss2015-2019-MI-descriptive"
r; t=0.00 9:51:17

.         local date "`year'-`month'-`day'"
r; t=0.00 9:51:17

.         local tag "`pgm'.do by `who' on `date'"
r; t=0.00 9:51:17

.         display "The tag is: `tag'"
The tag is: 20250823-brfss2015-2019-MI-descriptive.do by Ben Cheung on 2025-08-2
> 7
r; t=0.00 9:51:17

. 
. 
.         local today : di %tdCYND daily("$S_DATE", "DMY")
r; t=0.00 9:51:17

. 
.         *Check for installed packages
.         capture which scheme-plottig.scheme
r; t=0.00 9:51:17

.         if _rc!=0 ssc install blidschemes
r; t=0.00 9:51:17

. 
.         capture which confirmdir
r; t=0.00 9:51:17

.         if _rc!=0 ssc install confirmdir
r; t=0.00 9:51:17

. 
.         capture which estout
r; t=0.00 9:51:17

.         if _rc!=0 ssc install estout
r; t=0.00 9:51:17

. 
.         capture which sum2docx
r; t=0.00 9:51:17

.         if _rc!=0 ssc install sum2docx
r; t=0.00 9:51:17

. 
.         capture which tab2xl
r; t=0.00 9:51:17

.         if _rc!=0 net install http://www.stata.com/users/kcrow/tab2xl, replace
r; t=0.00 9:51:17

. 
.         set scheme plottig
r; t=0.00 9:51:17

. 
. 
. 
. 
. local folder "Data Ready for Analysis"
r; t=0.00 9:51:17

. confirmdir "`folder'"
r; t=0.00 9:51:17

. if `r(confirmdir)'==170 {
.         mkdir "`folder'"
r; t=0.00 9:51:17
. }
r; t=0.00 9:51:17

. 
. local folder "graph"
r; t=0.00 9:51:17

. confirmdir "`folder'"
r; t=0.00 9:51:17

. if `r(confirmdir)'==170 {
.         mkdir "`folder'"
r; t=0.00 9:51:17
. }
r; t=0.00 9:51:17

. 
. local folder "model"
r; t=0.00 9:51:17

. confirmdir "`folder'"
r; t=0.00 9:51:17

. if `r(confirmdir)'==170 {
.         mkdir "`folder'"
r; t=0.00 9:51:17
. }
r; t=0.00 9:51:17

. 
. local folder "summary"
r; t=0.00 9:51:17

. confirmdir "`folder'"
r; t=0.00 9:51:17

. if `r(confirmdir)'==170 {
.         mkdir "`folder'"
r; t=0.00 9:51:17
. }
r; t=0.00 9:51:17

. 
. 
. program define pe
  1.         version 15.1
  2.         if `"`0'"' != "" {
  3.                 display as input `"`0'"'
  4.                 `0'
  5.                 display("")
  6.         }
  7. end
r; t=0.00 9:51:17

. 
. 
. 
. // #1. Data management
. use "Data Ready for Analysis/brfss2015-MI-pathanalysis-sem-flong.dta",clear //
> Note that I use the multiply imputed dataset without dropping missing cases of
>  the two health outcomes
r; t=0.31 9:51:18

. 
. mi append using "Data Ready for Analysis/brfss2017-MI-pathanalysis-sem-flong.d
> ta"
r; t=61.24 9:52:19

. 
. mi append using "Data Ready for Analysis/brfss2019-MI-pathanalysis-sem-flong.d
> ta"
r; t=98.66 9:53:58

. 
. mi passive: gen yearsold2=yearsold*yearsold
m=0:
m=1:
m=2:
m=3:
m=4:
m=5:
m=6:
m=7:
m=8:
m=9:
m=10:
r; t=84.67 9:55:22

. 
. 
. 
. 
. 
. 
. 
. 
. mi estimate: svy, subpop(if male==0): proportion exe_aerobic_advice

Multiple-imputation estimates   Imputations     =           10
Survey: Proportion estimation   Number of obs   =    1,219,245

Number of strata  =     1,726   Population size =  657,058,941
Number of PSUs    = 1,219,245   Subpop. no. obs =      689,620
                                Subpop. size    =  339,643,007
                                Average RVI     =       2.2248
                                Largest FMI     =       0.8172
                                Complete DF     =      1217519
DF adjustment:   Small sample   DF:     min     =        14.30
                                        avg     =        15.40
Within VCE type:   Linearized           max     =        16.83

      _prop_1: exe_aerobic_advice = 0NoAerobicExercise
      _prop_2: exe_aerobic_advice = 1Under149MinEx
      _prop_3: exe_aerobic_advice = 2150+MinEx

--------------------------------------------------------------------
                   |                                   Normal
                   | Proportion   Std. Err.     [95% Conf. Interval]
-------------------+------------------------------------------------
           _prop_1 |   .3016729   .0022933      .2968307    .3065151
           _prop_2 |   .1994901   .0023125      .1945399    .2044403
           _prop_3 |    .498837   .0027227      .4930356    .5046384
--------------------------------------------------------------------
Note: 1 stratum omitted because it contains no subpopulation members.
r; t=106.28 9:57:09

. 
. mi estimate: svy, subpop(if male==0): proportion exe_strength_advice

Multiple-imputation estimates   Imputations     =           10
Survey: Proportion estimation   Number of obs   =    1,219,245

Number of strata  =     1,726   Population size =  657,058,941
Number of PSUs    = 1,219,245   Subpop. no. obs =      689,620
                                Subpop. size    =  339,643,007
                                Average RVI     =       0.4349
                                Largest FMI     =       0.3169
                                Complete DF     =      1217519
DF adjustment:   Small sample   DF:     min     =        97.96
                                        avg     =        97.96
Within VCE type:   Linearized           max     =        97.96

      _prop_1: exe_strength_advice = 0No
      _prop_2: exe_strength_advice = 1Yes

---------------------------------------------------------------------
                    |                                   Normal
                    | Proportion   Std. Err.     [95% Conf. Interval]
--------------------+------------------------------------------------
            _prop_1 |   .7378758   .0013725      .7351521    .7405996
            _prop_2 |   .2621242   .0013725      .2594004    .2648479
---------------------------------------------------------------------
Note: 1 stratum omitted because it contains no subpopulation members.
r; t=100.37 9:58:49

. 
. mi estimate: svy, subpop(if male==0): proportion praceethnicgp_multi_sim

Multiple-imputation estimates   Imputations     =           10
Survey: Proportion estimation   Number of obs   =    1,219,245

Number of strata  =     1,726   Population size =  657,058,941
Number of PSUs    = 1,219,245   Subpop. no. obs =      689,620
                                Subpop. size    =  339,643,007
                                Average RVI     =       0.0330
                                Largest FMI     =       0.0604
                                Complete DF     =      1217519
DF adjustment:   Small sample   DF:     min     =     2,524.95
                                        avg     =    12,093.77
Within VCE type:   Linearized           max     =    24,081.05

      _prop_1: praceethnicgp_multi_sim = 0nHWhite_Single
      _prop_2: praceethnicgp_multi_sim = 1nHwhite_Multi
      _prop_3: praceethnicgp_multi_sim = 2NHPI_Single
      _prop_4: praceethnicgp_multi_sim = 3NHPI_Multi
      _prop_5: praceethnicgp_multi_sim = 4Other

-------------------------------------------------------------------------
                        |                                   Normal
                        | Proportion   Std. Err.     [95% Conf. Interval]
------------------------+------------------------------------------------
                _prop_1 |   .6575545   .0013658      .6548775    .6602315
                _prop_2 |   .0052947   .0001456      .0050093    .0055802
                _prop_3 |   .0028239   .0001401      .0025493    .0030985
                _prop_4 |   .0007265   .0000443      .0006398    .0008133
                _prop_5 |   .3336003   .0013693      .3309163    .3362842
-------------------------------------------------------------------------
Note: 1 stratum omitted because it contains no subpopulation members.
r; t=123.95 10:00:53

. 
. mi estimate: svy, subpop(if male==0): mean yearsold

Multiple-imputation estimates   Imputations     =           10
Survey: Mean estimation         Number of obs   =    1,219,245

Number of strata  =     1,726   Population size =  657,058,941
Number of PSUs    = 1,219,245   Subpop. no. obs =      689,620
                                Subpop. size    =  339,643,007
                                Average RVI     =       0.0000
                                Largest FMI     =       0.0000
                                Complete DF     =      1217519
DF adjustment:   Small sample   DF:     min     = 1,217,517.00
                                        avg     = 1,217,517.00
Within VCE type:   Linearized           max     = 1,217,517.00

--------------------------------------------------------------
             |       Mean   Std. Err.     [95% Conf. Interval]
-------------+------------------------------------------------
    yearsold |   51.90551   .0406251      51.82589    51.98514
--------------------------------------------------------------
Note: 1 stratum omitted because it contains no subpopulation
      members.
r; t=90.14 10:02:23

. 
. *Reference: https://stats.stackexchange.com/questions/120097/is-it-possible-to
> -manually-calculate-standard-deviation-for-a-multiply-imputed-s
. mi estimate (sd:sqrt(_b[yearsold2] - _b[yearsold]*_b[yearsold])): svy, subpop(
> if male==0): mean yearsold yearsold2

Multiple-imputation estimates   Imputations     =           10
Survey: Mean estimation         Number of obs   =    1,219,245

Number of strata  =     1,726   Population size =  657,058,941
Number of PSUs    = 1,219,245   Subpop. no. obs =      689,620
                                Subpop. size    =  339,643,007
                                Average RVI     =       0.0000
                                Largest FMI     =       0.0000
                                Complete DF     =      1217519
DF adjustment:   Small sample   DF:     min     = 1,217,517.00
                                        avg     = 1,217,517.00
Within VCE type:   Linearized           max     = 1,217,517.00

--------------------------------------------------------------
             |       Mean   Std. Err.     [95% Conf. Interval]
-------------+------------------------------------------------
    yearsold |   51.90551   .0406251      51.82589    51.98514
   yearsold2 |   2954.634   4.216344       2946.37    2962.898
--------------------------------------------------------------

Transformations                 Average RVI     =       0.0000
                                Largest FMI     =       0.0000
                                Complete DF     =      1217519
DF adjustment:   Small sample   DF:     min     = 1,217,517.00
                                        avg     = 1,217,517.00
Within VCE type:   Linearized           max     = 1,217,517.00

           sd: sqrt(_b[yearsold2] - _b[yearsold]*_b[yearsold])

--------------------------------------------------------------
             |       Mean   Std. Err.     [95% Conf. Interval]
-------------+------------------------------------------------
          sd |   16.13852   .0192139      16.10086    16.17618
--------------------------------------------------------------
Note: 1 stratum omitted because it contains no subpopulation
      members.
r; t=88.39 10:03:51

. mi testtransform sd
note: assuming equal fractions of missing information

           sd: sqrt(_b[yearsold2] - _b[yearsold]*_b[yearsold])

 ( 1)  sd = 0

       F(  1, 1.2e+06) = 7.1e+05
            Prob > F =    0.0000
r; t=0.07 10:03:51

. 
. mi estimate: svy, subpop(if male==0): proportion educ

Multiple-imputation estimates   Imputations     =           10
Survey: Proportion estimation   Number of obs   =    1,219,245

Number of strata  =     1,726   Population size =  657,058,941
Number of PSUs    = 1,219,245   Subpop. no. obs =      689,620
                                Subpop. size    =  339,643,007
                                Average RVI     =       0.0066
                                Largest FMI     =       0.0164
                                Complete DF     =      1217519
DF adjustment:   Small sample   DF:     min     =    32,639.25
                                        avg     =   203,598.40
Within VCE type:   Linearized           max     =   361,277.58

      _prop_1: educ = 0BelowHighSch
      _prop_2: educ = 1HighSch
      _prop_3: educ = 2SomeCollege
      _prop_4: educ = 3CollegeorAbove

--------------------------------------------------------------
             |                                   Normal
             | Proportion   Std. Err.     [95% Conf. Interval]
-------------+------------------------------------------------
     _prop_1 |   .1320796   .0010559      .1300101    .1341491
     _prop_2 |   .2567458   .0011209      .2545487    .2589428
     _prop_3 |   .3125565   .0012258       .310154    .3149591
     _prop_4 |   .2986181   .0010994      .2964633    .3007728
--------------------------------------------------------------
Note: 1 stratum omitted because it contains no subpopulation
      members.
r; t=117.02 10:05:48

. 
. mi estimate: svy, subpop(if male==0): proportion hhincome

Multiple-imputation estimates   Imputations     =           10
Survey: Proportion estimation   Number of obs   =    1,219,245

Number of strata  =     1,726   Population size =  657,058,941
Number of PSUs    = 1,219,245   Subpop. no. obs =      689,620
                                Subpop. size    =  339,643,007
                                Average RVI     =       0.1593
                                Largest FMI     =       0.2213
                                Complete DF     =      1217519
DF adjustment:   Small sample   DF:     min     =       197.55
                                        avg     =       684.93
Within VCE type:   Linearized           max     =     1,914.89

      _prop_1: hhincome = 0_<10000
      _prop_2: hhincome = 1_10000-14999
      _prop_3: hhincome = 2_15000-19999
      _prop_4: hhincome = 3_20000-24999
      _prop_5: hhincome = 4_25000-34999
      _prop_6: hhincome = 5_35000-49999
      _prop_7: hhincome = 6_50000-74999
      _prop_8: hhincome = 7_75000+

--------------------------------------------------------------
             |                                   Normal
             | Proportion   Std. Err.     [95% Conf. Interval]
-------------+------------------------------------------------
     _prop_1 |   .0646281   .0007716      .0631091    .0661471
     _prop_2 |   .0592074   .0006624      .0579082    .0605065
     _prop_3 |    .083698   .0007559      .0822154    .0851805
     _prop_4 |   .0973139   .0008516      .0956389     .098989
     _prop_5 |   .1085365   .0009286      .1067058    .1103672
     _prop_6 |   .1288532   .0009429      .1269937    .1307127
     _prop_7 |    .143997   .0009618      .1421064    .1458876
     _prop_8 |   .3137659   .0013027      .3112032    .3163286
--------------------------------------------------------------
Note: 1 stratum omitted because it contains no subpopulation
      members.
r; t=117.81 10:07:46

. 
. mi estimate: svy, subpop(if male==0): proportion censusregion_sim

Multiple-imputation estimates   Imputations     =           10
Survey: Proportion estimation   Number of obs   =    1,219,245

Number of strata  =     1,726   Population size =  657,058,941
Number of PSUs    = 1,219,245   Subpop. no. obs =      689,620
                                Subpop. size    =  339,643,007
                                Average RVI     =       0.0000
                                Largest FMI     =       0.0000
                                Complete DF     =      1217519
DF adjustment:   Small sample   DF:     min     = 1,217,517.00
                                        avg     = 1,217,517.00
Within VCE type:   Linearized           max     = 1,217,517.00

      _prop_1: censusregion_sim = 0Eastern
      _prop_2: censusregion_sim = 1West
      _prop_3: censusregion_sim = 2Hawaii
      _prop_4: censusregion_sim = 3Guam

------------------------------------------------------------------
                 |                                   Normal
                 | Proportion   Std. Err.     [95% Conf. Interval]
-----------------+------------------------------------------------
         _prop_1 |   .7649712   .0010094      .7629929    .7669496
         _prop_2 |   .2302231   .0010104      .2282427    .2322034
         _prop_3 |   .0044157   .0000497      .0043183    .0045131
         _prop_4 |     .00039   .0000102        .00037    .0004101
------------------------------------------------------------------
Note: 1 stratum omitted because it contains no subpopulation
      members.
r; t=100.96 10:09:27

. 
. 
. 
. 
. 
. 
. 
. mi estimate: svy, subpop(if male==1): proportion exe_aerobic_advice

Multiple-imputation estimates   Imputations     =           10
Survey: Proportion estimation   Number of obs   =    1,219,246

Number of strata  =     1,727   Population size =  657,061,651
Number of PSUs    = 1,219,246   Subpop. no. obs =      529,626
                                Subpop. size    =  317,418,644
                                Average RVI     =       2.0418
                                Largest FMI     =       0.8222
                                Complete DF     =      1217519
DF adjustment:   Small sample   DF:     min     =        14.11
                                        avg     =        16.16
Within VCE type:   Linearized           max     =        18.97

      _prop_1: exe_aerobic_advice = 0NoAerobicExercise
      _prop_2: exe_aerobic_advice = 1Under149MinEx
      _prop_3: exe_aerobic_advice = 2150+MinEx

--------------------------------------------------------------------
                   |                                   Normal
                   | Proportion   Std. Err.     [95% Conf. Interval]
-------------------+------------------------------------------------
           _prop_1 |   .2832669   .0022926       .278468    .2880659
           _prop_2 |   .1958475   .0023219      .1909095    .2007856
           _prop_3 |   .5208855   .0031212      .5141961     .527575
--------------------------------------------------------------------
Note: Strata with single sampling unit centered at overall mean.
r; t=98.78 10:11:06

. 
. mi estimate: svy, subpop(if male==1): proportion exe_strength_advice

Multiple-imputation estimates   Imputations     =           10
Survey: Proportion estimation   Number of obs   =    1,219,246

Number of strata  =     1,727   Population size =  657,061,651
Number of PSUs    = 1,219,246   Subpop. no. obs =      529,626
                                Subpop. size    =  317,418,644
                                Average RVI     =       0.5910
                                Largest FMI     =       0.3899
                                Complete DF     =      1217519
DF adjustment:   Small sample   DF:     min     =        65.22
                                        avg     =        65.22
Within VCE type:   Linearized           max     =        65.22

      _prop_1: exe_strength_advice = 0No
      _prop_2: exe_strength_advice = 1Yes

---------------------------------------------------------------------
                    |                                   Normal
                    | Proportion   Std. Err.     [95% Conf. Interval]
--------------------+------------------------------------------------
            _prop_1 |   .6730681   .0016621      .6697489    .6763873
            _prop_2 |   .3269319   .0016621      .3236127    .3302511
---------------------------------------------------------------------
Note: Strata with single sampling unit centered at overall mean.
r; t=96.58 10:12:43

. 
. mi estimate: svy, subpop(if male==1): proportion praceethnicgp_multi_sim

Multiple-imputation estimates   Imputations     =           10
Survey: Proportion estimation   Number of obs   =    1,219,246

Number of strata  =     1,727   Population size =  657,061,651
Number of PSUs    = 1,219,246   Subpop. no. obs =      529,626
                                Subpop. size    =  317,418,644
                                Average RVI     =       0.0383
                                Largest FMI     =       0.0620
                                Complete DF     =      1217519
DF adjustment:   Small sample   DF:     min     =     2,396.43
                                        avg     =     5,613.67
Within VCE type:   Linearized           max     =    14,386.00

      _prop_1: praceethnicgp_multi_sim = 0nHWhite_Single
      _prop_2: praceethnicgp_multi_sim = 1nHwhite_Multi
      _prop_3: praceethnicgp_multi_sim = 2NHPI_Single
      _prop_4: praceethnicgp_multi_sim = 3NHPI_Multi
      _prop_5: praceethnicgp_multi_sim = 4Other

-------------------------------------------------------------------------
                        |                                   Normal
                        | Proportion   Std. Err.     [95% Conf. Interval]
------------------------+------------------------------------------------
                _prop_1 |   .6579356   .0014945      .6550051    .6608662
                _prop_2 |   .0062576   .0001719      .0059207    .0065946
                _prop_3 |   .0035581    .000162      .0032405    .0038757
                _prop_4 |   .0007072   .0000453      .0006185    .0007959
                _prop_5 |   .3315415   .0014969      .3286064    .3344765
-------------------------------------------------------------------------
Note: Strata with single sampling unit centered at overall mean.
r; t=116.99 10:14:40

. 
. mi estimate: svy, subpop(if male==1): mean yearsold

Multiple-imputation estimates   Imputations     =           10
Survey: Mean estimation         Number of obs   =    1,219,246

Number of strata  =     1,727   Population size =  657,061,651
Number of PSUs    = 1,219,246   Subpop. no. obs =      529,626
                                Subpop. size    =  317,418,644
                                Average RVI     =       0.0000
                                Largest FMI     =       0.0000
                                Complete DF     =      1217519
DF adjustment:   Small sample   DF:     min     = 1,217,517.00
                                        avg     = 1,217,517.00
Within VCE type:   Linearized           max     = 1,217,517.00

--------------------------------------------------------------
             |       Mean   Std. Err.     [95% Conf. Interval]
-------------+------------------------------------------------
    yearsold |    50.3314   .0424845      50.24813    50.41466
--------------------------------------------------------------
Note: Strata with single sampling unit centered at overall
      mean.
r; t=77.30 10:15:57

. 
. *Reference: https://stats.stackexchange.com/questions/120097/is-it-possible-to
> -manually-calculate-standard-deviation-for-a-multiply-imputed-s
. mi estimate (sd:sqrt(_b[yearsold2] - _b[yearsold]*_b[yearsold])): svy, subpop(
> if male==1): mean yearsold yearsold2

Multiple-imputation estimates   Imputations     =           10
Survey: Mean estimation         Number of obs   =    1,219,246

Number of strata  =     1,727   Population size =  657,061,651
Number of PSUs    = 1,219,246   Subpop. no. obs =      529,626
                                Subpop. size    =  317,418,644
                                Average RVI     =       0.0000
                                Largest FMI     =       0.0000
                                Complete DF     =      1217519
DF adjustment:   Small sample   DF:     min     = 1,217,517.00
                                        avg     = 1,217,517.00
Within VCE type:   Linearized           max     = 1,217,517.00

--------------------------------------------------------------
             |       Mean   Std. Err.     [95% Conf. Interval]
-------------+------------------------------------------------
    yearsold |    50.3314   .0424845      50.24813    50.41466
   yearsold2 |   2780.041   4.331955      2771.551    2788.532
--------------------------------------------------------------

Transformations                 Average RVI     =       0.0000
                                Largest FMI     =       0.0000
                                Complete DF     =      1217519
DF adjustment:   Small sample   DF:     min     = 1,217,517.00
                                        avg     = 1,217,517.00
Within VCE type:   Linearized           max     = 1,217,517.00

           sd: sqrt(_b[yearsold2] - _b[yearsold]*_b[yearsold])

--------------------------------------------------------------
             |       Mean   Std. Err.     [95% Conf. Interval]
-------------+------------------------------------------------
          sd |   15.70961    .020502      15.66943    15.74979
--------------------------------------------------------------
Note: Strata with single sampling unit centered at overall
      mean.
r; t=80.14 10:17:17

. mi testtransform sd
note: assuming equal fractions of missing information

           sd: sqrt(_b[yearsold2] - _b[yearsold]*_b[yearsold])

 ( 1)  sd = 0

       F(  1, 1.2e+06) = 5.9e+05
            Prob > F =    0.0000
r; t=0.00 10:17:17

. 
. mi estimate: svy, subpop(if male==1): proportion educ

Multiple-imputation estimates   Imputations     =           10
Survey: Proportion estimation   Number of obs   =    1,219,246

Number of strata  =     1,727   Population size =  657,061,651
Number of PSUs    = 1,219,246   Subpop. no. obs =      529,626
                                Subpop. size    =  317,418,644
                                Average RVI     =       0.0030
                                Largest FMI     =       0.0047
                                Complete DF     =      1217519
DF adjustment:   Small sample   DF:     min     =   301,062.34
                                        avg     =   394,824.83
Within VCE type:   Linearized           max     =   441,720.12

      _prop_1: educ = 0BelowHighSch
      _prop_2: educ = 1HighSch
      _prop_3: educ = 2SomeCollege
      _prop_4: educ = 3CollegeorAbove

--------------------------------------------------------------
             |                                   Normal
             | Proportion   Std. Err.     [95% Conf. Interval]
-------------+------------------------------------------------
     _prop_1 |    .144116   .0011959      .1417721    .1464599
     _prop_2 |   .2754427     .00124      .2730122    .2778731
     _prop_3 |   .2848562    .001295      .2823182    .2873943
     _prop_4 |   .2955851    .001148      .2933351    .2978352
--------------------------------------------------------------
Note: Strata with single sampling unit centered at overall
      mean.
r; t=102.59 10:19:00

. 
. mi estimate: svy, subpop(if male==1): proportion hhincome

Multiple-imputation estimates   Imputations     =           10
Survey: Proportion estimation   Number of obs   =    1,219,246

Number of strata  =     1,727   Population size =  657,061,651
Number of PSUs    = 1,219,246   Subpop. no. obs =      529,626
                                Subpop. size    =  317,418,644
                                Average RVI     =       0.1415
                                Largest FMI     =       0.2045
                                Complete DF     =      1217519
DF adjustment:   Small sample   DF:     min     =       230.48
                                        avg     =       542.65
Within VCE type:   Linearized           max     =       920.15

      _prop_1: hhincome = 0_<10000
      _prop_2: hhincome = 1_10000-14999
      _prop_3: hhincome = 2_15000-19999
      _prop_4: hhincome = 3_20000-24999
      _prop_5: hhincome = 4_25000-34999
      _prop_6: hhincome = 5_35000-49999
      _prop_7: hhincome = 6_50000-74999
      _prop_8: hhincome = 7_75000+

--------------------------------------------------------------
             |                                   Normal
             | Proportion   Std. Err.     [95% Conf. Interval]
-------------+------------------------------------------------
     _prop_1 |   .0438262   .0006712      .0425062    .0451462
     _prop_2 |   .0423386   .0006104      .0411406    .0435365
     _prop_3 |   .0648114   .0007642      .0633107     .066312
     _prop_4 |   .0826984   .0008682      .0809934    .0844034
     _prop_5 |   .0959691   .0009382      .0941205    .0978176
     _prop_6 |    .134108   .0010346       .132073    .1361429
     _prop_7 |   .1559936   .0010532      .1539245    .1580627
     _prop_8 |   .3802547   .0014278      .3774516    .3830579
--------------------------------------------------------------
Note: Strata with single sampling unit centered at overall
      mean.
r; t=115.64 10:20:55

. 
. mi estimate: svy, subpop(if male==1): proportion censusregion_sim

Multiple-imputation estimates   Imputations     =           10
Survey: Proportion estimation   Number of obs   =    1,219,246

Number of strata  =     1,727   Population size =  657,061,651
Number of PSUs    = 1,219,246   Subpop. no. obs =      529,626
                                Subpop. size    =  317,418,644
                                Average RVI     =       0.0000
                                Largest FMI     =       0.0000
                                Complete DF     =      1217519
DF adjustment:   Small sample   DF:     min     = 1,217,517.00
                                        avg     = 1,217,517.00
Within VCE type:   Linearized           max     = 1,217,517.00

      _prop_1: censusregion_sim = 0Eastern
      _prop_2: censusregion_sim = 1West
      _prop_3: censusregion_sim = 2Hawaii
      _prop_4: censusregion_sim = 3Guam

------------------------------------------------------------------
                 |                                   Normal
                 | Proportion   Std. Err.     [95% Conf. Interval]
-----------------+------------------------------------------------
         _prop_1 |   .7573945   .0010406      .7553549    .7594341
         _prop_2 |   .2375254   .0010405      .2354861    .2395647
         _prop_3 |   .0046538   .0000536      .0045487    .0047589
         _prop_4 |   .0004262    .000011      .0004046    .0004478
------------------------------------------------------------------
Note: Strata with single sampling unit centered at overall mean.
r; t=100.07 10:22:35

. 
. 
. 
. 
. 
. 
. 
. mi estimate: svy: proportion exe_aerobic_advice

Multiple-imputation estimates   Imputations     =           10
Survey: Proportion estimation   Number of obs   =    1,219,246

Number of strata  =     1,727   Population size =  657,061,651
Number of PSUs    = 1,219,246
                                Average RVI     =       6.1660
                                Largest FMI     =       0.8944
                                Complete DF     =      1217519
DF adjustment:   Small sample   DF:     min     =        11.68
                                        avg     =        12.23
Within VCE type:   Linearized           max     =        13.11

      _prop_1: exe_aerobic_advice = 0NoAerobicExercise
      _prop_2: exe_aerobic_advice = 1Under149MinEx
      _prop_3: exe_aerobic_advice = 2150+MinEx

--------------------------------------------------------------------
                   |                                   Normal
                   | Proportion   Std. Err.     [95% Conf. Interval]
-------------------+------------------------------------------------
           _prop_1 |   .2927812   .0021011      .2882459    .2973165
           _prop_2 |   .1977304   .0021344       .193075    .2023859
           _prop_3 |   .5094884   .0027217      .5035403    .5154364
--------------------------------------------------------------------
Note: Strata with single sampling unit centered at overall mean.
r; t=98.20 10:24:14

. 
. mi estimate: svy: proportion exe_strength_advice

Multiple-imputation estimates   Imputations     =           10
Survey: Proportion estimation   Number of obs   =    1,219,246

Number of strata  =     1,727   Population size =  657,061,651
Number of PSUs    = 1,219,246
                                Average RVI     =       0.7563
                                Largest FMI     =       0.4527
                                Complete DF     =      1217519
DF adjustment:   Small sample   DF:     min     =        48.53
                                        avg     =        48.53
Within VCE type:   Linearized           max     =        48.53

      _prop_1: exe_strength_advice = 0No
      _prop_2: exe_strength_advice = 1Yes

---------------------------------------------------------------------
                    |                                   Normal
                    | Proportion   Std. Err.     [95% Conf. Interval]
--------------------+------------------------------------------------
            _prop_1 |    .706568   .0011545      .7042473    .7088887
            _prop_2 |    .293432   .0011545      .2911113    .2957527
---------------------------------------------------------------------
Note: Strata with single sampling unit centered at overall mean.
r; t=91.24 10:25:45

. 
. mi estimate: svy: proportion praceethnicgp_multi_sim

Multiple-imputation estimates   Imputations     =           10
Survey: Proportion estimation   Number of obs   =    1,219,246

Number of strata  =     1,727   Population size =  657,061,651
Number of PSUs    = 1,219,246
                                Average RVI     =       0.0326
                                Largest FMI     =       0.0406
                                Complete DF     =      1217519
DF adjustment:   Small sample   DF:     min     =     5,521.93
                                        avg     =    10,485.42
Within VCE type:   Linearized           max     =    14,606.62

      _prop_1: praceethnicgp_multi_sim = 0nHWhite_Single
      _prop_2: praceethnicgp_multi_sim = 1nHwhite_Multi
      _prop_3: praceethnicgp_multi_sim = 2NHPI_Single
      _prop_4: praceethnicgp_multi_sim = 3NHPI_Multi
      _prop_5: praceethnicgp_multi_sim = 4Other

-------------------------------------------------------------------------
                        |                                   Normal
                        | Proportion   Std. Err.     [95% Conf. Interval]
------------------------+------------------------------------------------
                _prop_1 |   .6577386   .0009764      .6558248    .6596524
                _prop_2 |   .0057599    .000111      .0055422    .0059776
                _prop_3 |   .0031786   .0001059       .002971    .0033862
                _prop_4 |   .0007172   .0000315      .0006554    .0007789
                _prop_5 |   .3326057   .0009802      .3306844     .334527
-------------------------------------------------------------------------
Note: Strata with single sampling unit centered at overall mean.
r; t=161.59 10:28:26

. 
. mi estimate: svy: mean yearsold

Multiple-imputation estimates   Imputations     =           10
Survey: Mean estimation         Number of obs   =    1,219,246

Number of strata  =     1,727   Population size =  657,061,651
Number of PSUs    = 1,219,246
                                Average RVI     =       0.0000
                                Largest FMI     =       0.0000
                                Complete DF     =      1217519
DF adjustment:   Small sample   DF:     min     = 1,217,517.00
                                        avg     = 1,217,517.00
Within VCE type:   Linearized           max     = 1,217,517.00

--------------------------------------------------------------
             |       Mean   Std. Err.     [95% Conf. Interval]
-------------+------------------------------------------------
    yearsold |   51.14508   .0288631      51.08851    51.20165
--------------------------------------------------------------
Note: Strata with single sampling unit centered at overall
      mean.
r; t=121.53 10:30:28

. 
. *Reference: https://stats.stackexchange.com/questions/120097/is-it-possible-to
> -manually-calculate-standard-deviation-for-a-multiply-imputed-s
. mi estimate (sd:sqrt(_b[yearsold2] - _b[yearsold]*_b[yearsold])): svy: mean ye
> arsold yearsold2

Multiple-imputation estimates   Imputations     =           10
Survey: Mean estimation         Number of obs   =    1,219,246

Number of strata  =     1,727   Population size =  657,061,651
Number of PSUs    = 1,219,246
                                Average RVI     =       0.0000
                                Largest FMI     =       0.0000
                                Complete DF     =      1217519
DF adjustment:   Small sample   DF:     min     = 1,217,517.00
                                        avg     = 1,217,517.00
Within VCE type:   Linearized           max     = 1,217,517.00

--------------------------------------------------------------
             |       Mean   Std. Err.     [95% Conf. Interval]
-------------+------------------------------------------------
    yearsold |   51.14508   .0288631      51.08851    51.20165
   yearsold2 |    2870.29   2.965995      2864.477    2876.104
--------------------------------------------------------------

Transformations                 Average RVI     =       0.0000
                                Largest FMI     =       0.0000
                                Complete DF     =      1217519
DF adjustment:   Small sample   DF:     min     = 1,217,517.00
                                        avg     = 1,217,517.00
Within VCE type:   Linearized           max     = 1,217,517.00

           sd: sqrt(_b[yearsold2] - _b[yearsold]*_b[yearsold])

--------------------------------------------------------------
             |       Mean   Std. Err.     [95% Conf. Interval]
-------------+------------------------------------------------
          sd |   15.95217   .0139875      15.92475    15.97958
--------------------------------------------------------------
Note: Strata with single sampling unit centered at overall
      mean.
r; t=121.37 10:32:29

. mi testtransform sd
note: assuming equal fractions of missing information

           sd: sqrt(_b[yearsold2] - _b[yearsold]*_b[yearsold])

 ( 1)  sd = 0

       F(  1, 1.2e+06) = 1.3e+06
            Prob > F =    0.0000
r; t=0.00 10:32:29

. 
. mi estimate: svy: proportion educ

Multiple-imputation estimates   Imputations     =           10
Survey: Proportion estimation   Number of obs   =    1,219,246

Number of strata  =     1,727   Population size =  657,061,651
Number of PSUs    = 1,219,246
                                Average RVI     =       0.0047
                                Largest FMI     =       0.0120
                                Complete DF     =      1217519
DF adjustment:   Small sample   DF:     min     =    59,585.12
                                        avg     =   286,327.72
Within VCE type:   Linearized           max     =   451,684.07

      _prop_1: educ = 0BelowHighSch
      _prop_2: educ = 1HighSch
      _prop_3: educ = 2SomeCollege
      _prop_4: educ = 3CollegeorAbove

--------------------------------------------------------------
             |                                   Normal
             | Proportion   Std. Err.     [95% Conf. Interval]
-------------+------------------------------------------------
     _prop_1 |   .1378942   .0007927      .1363406    .1394479
     _prop_2 |    .265778   .0008324      .2641465    .2674095
     _prop_3 |   .2991748   .0008904      .2974297      .30092
     _prop_4 |   .2971529   .0007911      .2956025    .2987033
--------------------------------------------------------------
Note: Strata with single sampling unit centered at overall
      mean.
r; t=106.61 10:34:18

. 
. mi estimate: svy: proportion hhincome

Multiple-imputation estimates   Imputations     =           10
Survey: Proportion estimation   Number of obs   =    1,219,246

Number of strata  =     1,727   Population size =  657,061,651
Number of PSUs    = 1,219,246
                                Average RVI     =       0.1598
                                Largest FMI     =       0.1993
                                Complete DF     =      1217519
DF adjustment:   Small sample   DF:     min     =       242.29
                                        avg     =       569.15
Within VCE type:   Linearized           max     =     1,190.48

      _prop_1: hhincome = 0_<10000
      _prop_2: hhincome = 1_10000-14999
      _prop_3: hhincome = 2_15000-19999
      _prop_4: hhincome = 3_20000-24999
      _prop_5: hhincome = 4_25000-34999
      _prop_6: hhincome = 5_35000-49999
      _prop_7: hhincome = 6_50000-74999
      _prop_8: hhincome = 7_75000+

--------------------------------------------------------------
             |                                   Normal
             | Proportion   Std. Err.     [95% Conf. Interval]
-------------+------------------------------------------------
     _prop_1 |    .054579   .0005043      .0535879      .05557
     _prop_2 |   .0510583   .0004578      .0501596    .0519569
     _prop_3 |   .0745741   .0005356      .0735231     .075625
     _prop_4 |   .0902533   .0006262      .0890199    .0914868
     _prop_5 |   .1024653    .000629        .10123    .1037006
     _prop_6 |   .1313917   .0006788      .1300578    .1327257
     _prop_7 |   .1497924   .0007257      .1483642    .1512206
     _prop_8 |   .3458859   .0009625      .3439946    .3477771
--------------------------------------------------------------
Note: Strata with single sampling unit centered at overall
      mean.
r; t=124.92 10:36:23

. 
. mi estimate: svy: proportion censusregion_sim

Multiple-imputation estimates   Imputations     =           10
Survey: Proportion estimation   Number of obs   =    1,219,246

Number of strata  =     1,727   Population size =  657,061,651
Number of PSUs    = 1,219,246
                                Average RVI     =       0.0000
                                Largest FMI     =       0.0000
                                Complete DF     =      1217519
DF adjustment:   Small sample   DF:     min     = 1,217,517.00
                                        avg     = 1,217,517.00
Within VCE type:   Linearized           max     = 1,217,517.00

      _prop_1: censusregion_sim = 0Eastern
      _prop_2: censusregion_sim = 1West
      _prop_3: censusregion_sim = 2Hawaii
      _prop_4: censusregion_sim = 3Guam

------------------------------------------------------------------
                 |                                   Normal
                 | Proportion   Std. Err.     [95% Conf. Interval]
-----------------+------------------------------------------------
         _prop_1 |    .761311   .0005515        .76023     .762392
         _prop_2 |   .2337507   .0005518      .2326692    .2348323
         _prop_3 |   .0045307    .000027      .0044778    .0045837
         _prop_4 |   .0004075   6.29e-06      .0003952    .0004198
------------------------------------------------------------------
Note: Strata with single sampling unit centered at overall mean.
r; t=109.07 10:38:12

. 
. 
. 
. 
. 
. log close
